arXiv: 1503.00192vl [math-ph] 28 Feb 2015 


A COMPACTNESS LEMMA AND ITS APPLICATION TO THE 
EXISTENCE OF MINIMIZERS FOR THE LIQUID DROP MODEL 


RUPERT L. FRANK AND ELLIOTT H. LIEB 

Abstract. The ancient Gamow liquid drop model of nuclear energies has had a 
renewed life as an interesting problem in the calculus of variations: Find a set P C M 3 
with given volume A that minimizes the sum of its surface area and its Coulomb self 
energy. A ball minimizes the former and maximizes the latter, but the conjecture is 
that a ball is always a minimizer - when there is a minimizer. Even the existence of 
minimizers for this interesting geometric problem has not been shown in general. We 
prove the existence of the absolute minimizer (over all A) of the energy divided by 
A (the binding energy per particle). A second result of our work is a general method 
for showing the existence of optimal sets in geometric minimization problems, which 
we call the ‘method of the missing mass’. A third point is the extension of the pulling 
back compactness lemma [L3] from W 1,p to BV. 


1. Introduction and main results 

In this paper we are interested in some aspects of Gamow’s famous 1930 ‘liquid 
drop model’ jG]. This simple, but successful model of atomic nuclei predicts, among 
other things, 

(a) the spherical shape of nuclei, 

(b) the non-existence of nuclei with atomic number larger than some critical value, 
and 

(c) the existence of a nucleus with minimal binding energy per particle. 

Rather surprisingly, this model has received little attention in mathematics, and only 
recently have proofs appeared which rigorously derive properties (a) (partially) and 
(b). One of our goals in this paper will be to prove property (c). Before we describe 
these results in detail, we recall the model. 

In the liquid drop model a nucleus O C M 3 is assumed to have constant density, 
which we may assume to be equal to one. Thus, |0| = A is the number of nucleons 
(protons and neutrons) in the nucleus. Mathematically, it is not necessary to assume 
that this number is an integer. The binding energy of a nucleus is given by 

£(fi) := Per Q + D(Q), 
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where Per denotes the surface area of provided its boundary is smooth (see (12.11) 
for a definition of Peril for an arbitrary measurable set), and 

0<n ):=Iff ^ 

denotes the Coulomb repulsion energy. We dropped a volume term — eA, where e is the 
energy of ‘nuclear matter’ per unit volume in the ‘thermodynamic limit’. (This energy 
is computed using only the short-range strong forces, without the Coulomb repulsion.) 
This term only contributes a constant in our case, and need not be exhibited explicitly 
for our considerations, but stability of a real nucleus requires its inclusion; the total 
energy must be negative. The surface term Peril accounts for corrections to the 
volume term coming from the surface of the nucleus and acts as a surface tension. 
The Coulomb term D(Q) describes the repulsion between the protons in the nucleus. 
The charge (which we set equal to one in the definition of -D(ll)) is an effective charge 
proportional to the ratio of the proton number to A. We also neglected asymmetry 
and pairing terms. 

The minimal binding energy of a nucleus with mass number A > 0 is given by 

E(A) := inf{5(H) : |H| = A} . 

The following is rigorously known about the properties (a) and (b) mentioned before. 

(a) There is an A Cl > 0 such that E(A) is minimized by a ball for all 0 < A < A Cl 
and such that a ball does not minimize E(A) for A > A C1 . It is even known 
that balls are the unique minimizers for A < A C1 . 

(b) There is an A C2 > 0 such that E(A) has no minimizer for A > A C2 . 

Fact (a) was proved by Kniipfer and Muratov [ KM1 Thm. 3.2] (balls are minimizers 
for A < A Cl ] see also CD and by Bonacini and Cristoferi [ BC. Thm. 2.10] (balls are 
unique minimizers for A < A Cl and are not minimizers for A > A Cl ). Fact (b) is 
proved by Lu and Otto TO, Thm. 2] and their 2014 presentation of this work at the 
Fields Institute was partial motivation for our interest in the problem. 

Fact (b) has the physical interpretation that nuclei ‘fission’ spontaneously when 
they are too large. When the mass is slightly below the critical value A C2 they are 
still unstable against being impacted by a neutron - with the release of some notable 
amount of energy. In this model the reduction in energy is electrostatic (Coulomb), 
not nuclear, so one might properly refer to a ‘Coulomb bomb’ rather than to a ‘nuclear 
bomb’. 

It is conjectured ra, but not known, that A Cl = A C2 . 

In this paper we show that (see Theorem 13.21) : 

(c) There is a set hi* C R 3 such that 

E(A) 

A ~ |fi*| 


for any A > 0. 
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Physically, the set corresponds to an iron (Fe) nucleus which is the element with 
the greatest binding energy per particle. Our proof of (c) makes use of the Lu-Otto 
analysis in |LOj . From (c) one can easily deduce that (see Corollary 13.311 

A—>-oo y± |£z*| 

Recall from (b) that 

A C2 := sup {A > 0 : E(A) has a minimizer} < oo . 

We consider the question whether there is a minimizer for E(A C2 ). By definition there 
is a sequence of sets Q n with A n := |h2 n | —y A C2 which minimize E(A n ). The question 
is how these sets behave as n —> oo. ff they have a limit, this limit is a good candidate 
for a minimizer for E(A C2 ). On the other hand, the fact that there is no minimizer 
beyond A C2 means, in some sense, that the Coulomb energy becomes the dominant 
term at A C2 . In order to minimize the Coulomb energy, one might think that the sets 
Vt n become more and more elongated, which could cause the failure of a non-trivial 
limit. That this possible scenario does not occur is our next result (see Theorem 13.4p . 

(d) E(A C2 ) has a minimizer. 

Physically, a corresponding minimizing set Vt C2 corresponds to lead (Pb), a stable 
nucleus with the largest possible nucleon number. Our proof of (d) uses crucially 
bounds from [LO] for the proof of (b). 

We also prove the following (physically obvious) fact (see Theorem 13.111 . 

(e) If for a given A > 0 one has the strict binding inequality E(A) < E(A') + 
E(A — A') for all 0 < A' < A, then there is a minimizer for E(A). 

We emphasize that the strict binding inequality is only a sufficient condition for the 
existence of a minimizer. For instance in the situation of (d) it is probably violated. 

A second goal in this paper, besides the study of the liquid drop model, is to present 
a method to prove the existence of optimal sets in geometric minimization problems 
that involve the perimeter functional. We believe that our method is simpler and 
more direct than methods used, for instance, in [dPV] and has applications in related 
problems. 

The method that we are using may be called the ‘method of the missing mass ’ and 
originates in the solution of the HLS problem |L2j and also used in the paper [BNl 
Lemma 1.2 and Appendix], In this method each element of a minimizing sequence is 
decomposed into two terms, namely, a main piece with good convergence properties 
and a missing piece which vanishes in a weak sense. The key point is then to use this 
missing piece as a potential minimizer. Typically, this allows one to conclude that the 
limit of the main piece is a minimizer. This method is often, but not in this paper, 
combined with a strengthening of Fatou’s lemma which includes the missing piece; see 
|L2| and, for an extension to a wide class of convex functions, see (BLll . (The latter 
paper also contains a review of how this method is used in j|L2| and |BNj .) 
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The method of the missing mass works both in problems with and without trans¬ 
lation invariance. The minimization problems in |BN] is not translation invariant and 
the translation invariance in m is broken by symmetric decreasing rearrangement. 
In situations with translation invariance, however, as in the present paper, there is 
the additional difficulty of detecting a main piece which has a non-zero limit. This 
problem was overcome in |L3j : see (BL2l iFLL j for early applications of this pulling 
back compactness lemma to concrete minimization problems. More recent applications 
include |FLl IBFV] . 

As far as we know, the method of the missing mass has not been applied to geometric 
minimization problems in which the minimization is over sets. In this paper we identify 
a possible non-zero limit after translations in much the same way as in |L3j . where the 
case of the Sobolev space IF 1,p (M d ) with p > 1 was treated. Our case here corresponds 
to p — 1, where the space BV(W l ) of functions of bounded variation has to be used 
to get better compactness properties; see Proposition 12.11 What is considerably more 
involved than in the p = 2 case is to split the elements of the minimizing sequence 
into a main and a missing piece. This is accomplished in Lemma 12.21 and the method 
used has some similarities with the one used in (FLSSj in a different physical context. 

2. Compactness up to translations 

In this section we describe sequences of subsets of M. d with a uniform perimeter 
bound and we prove a pulling back compactness theorem, which extends that of [ L3l 
Sect. 3] to functions of bounded variation. 

If E n and E are measurable sets in W 1 , we say that E n —>■ E (globally) if \E n AE\ —>• 0 
and we say that E n —>■ E locally if for every compact set K c M. d , 

\(E n AE)DK\ -A 0. 

If E C is a measurable set, we denote by 

Per E = sup | J di v F dx : F 6 C' 0 1 (R (i ,R‘ i ), \F\ < lj (2.1) 

its perimeter in the sense of De Giorgi. Moreover, B r (a) denotes the open ball of 
radius r centered at a G M. d , and we abbreviate B r = £> r (0). 

The following proposition is our technical main result. As explained in the intro¬ 
duction, its proof relies on a technique from H. 

Proposition 2.1. Let ( E n ) be a sequence of measurable sets in W l with uniformly 
bounded perimeter. Then one of the following two alternatives occurs: 

(1) liny, ,^ |F/ n | 0 

(2) There is a set E with positive measure and a sequence (a*,) in W. d such that for 
a subsequence (nk) one has 

E nk — ak —» E locally. 

Moreover, \E\ < liniinf),^^ \E nk \ and Per E < liminf),^^Per E nk . 
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Proof. We begin by proving that for every u G W 1,1 (R <i ) and every r > 0, 

-1/d 

(|Vu| + CV' _1 |w|) dx > c [ sup \B r (a) fl {u ^ 0}| ) / \u\ dx 

\a£ R d 


( 2 . 2 ) 


for some constants c and C depending only on d. In fact, let x be a smooth function 
with support in B 1 and ||x||i = 1- We put Xa,r{x) = r~ d x({x — a)/r ) and compute 


|V(x a , r u)| dx < 


(|Xa,r||Vu| + \S7XaA\ U \) dx ■ 


On the other hand, by a Sobolev inquality and Holder inequality, 

(d-l)/d 

|V {Xa, r u)\ dx > C ( / \Xa, r u\ d/{d ~ 1) dx 


> c ( sup | B r (a) n{u^ 0}| 

\a€ R d 


-I/d 


| Xa,r II | dx . 


When we integrate both the upper and the lower bound on ||V(Xa,r-w)||i with respect 
to a we obtain (12.21) with C = ||Vx||i- 

By a standard argument [ AFP . Thm. 3.9] (12.21) extends to u G and there¬ 

fore, in particular, to u — Xe„, where (E n ) has uniformly bounded perimeter. For 
r = 1, we find that 


sup |-£>i(a) D E n \ > 

aSM d 


c| E ti 


P er E n + C | E n 


d 


If we assume that alternative (1) does not occur, we may pass to a subsequence and 
assume that inf n \E n \ > 0. Thus the right side of the previous inequality is bounded 
away from zero and so there are a n G such that 


inf \B i n (E n - a n )| = inf \Bi(a n ) D E n \ > 0 . 

n n 

We now apply the compactness result for sets of finite perimeter [AFP1 Thm. 3.39] to 
E n —a n . We infer that there is a set E of finite perimeter such that E n — a n ^-E locally. 
It is easy to see that local convergence implies that liminf \E n \ — liminf \E n — a n \ > 
| E\. Finally, lower semi-continuity of the perimeter, that is, liin inf Per E n = Per (E n — 
a n ) > Per A 1 , follows easily from the variational dehnition of the perimeter. □ 


Next, we describe sequences of sets which converge locally but not globally, because 
they lose measure in the limit. We show that such sets can be decomposed into a main 
piece, which does converge globally, and a remainder piece which disappears from any 
compact set. 


Lemma 2.2. Let (E n ) be a sequence of measurable sets in with finite perimeter such 
that E n —^ E locally for some set E C W 1 . Assume that 0 < |A| < liminf^oo |A„|. 
Then there is a sequence (r n ) in (0, oo) such that the sets 


F n := E n n B rn , 


G n := E n n ( R d \ B rn ) 
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satisfy 

and 

In particular, 


lim (Per E n — Per F n — Per G n ) = 0 

n—>oo 

F n —$■ E (globally) and G n —* 0 locally. 
lim \F n \ = \E\ and lim inf Per F n > Per E . 

n—> oo n—too 


Proof. Clearly, we have 

Per E n < Per F n + Per G n < Per E n + 2a rn ( dB Vn n E n ), 


(2.3) 

(2.4) 

(2.5) 


where a r denotes surface measure on the sphere dB r . Thus, for the bound (1 2. ‘ill on 
the perimeter, we have to prove that 


<T r n {dB rn C En) 0 (2.6) 

for a suitable choice of r n . 

In order to construct the r n we distinguish two cases. We hrst assume that E is 
(essentially) bounded, say, E C B R for some R > 0. Then E n —>• E locally implies 
that e n := | E n D (B 2R \ B R )| — > 0. We claim that for every n there is an r n E [/?, 2 R] 
with cr rn (dB rn D E n ) < £ n /R, proving (12.61) . In fact, if this was not the case, we had, 
by integration in spherical coordinates, 

r2R 

£ n < a r (dB r n E n ) dr = \E n n ( B 2 r \ B R ) \=£ n , 

Jr 

which is a contradiction. 

Now assume that E is (essentially) unbounded and define R n by \E n D B Rn \ = \E\. 
(This is well defined if \E n \ > \E\ , which we may assume after discarding finitely many 
En s.) We claim that Rn oo. In fact, if we had R nk —> R * for some R* < oo, then, 
using E n —)• E locally, 


1^1 = I E nk n B Rnk I = I E nk n B r , 1 1 + o(l) = \EnBR,\ + o(l) < 1^1 + o(l), 

which is a contradiction. Now choose R such that \E D B R \ = \E\/2. Then, since 
E n —)■ E locally, for all sufficiently large n, \E n D B R \ > \E\/4. After discarding 
finitely many n we may assume that R n > R. We claim that for every n there is an 
t n e [(R + R n )/2,R n \ with a rn (dB rn n E n ) < 3\E\/(2(R n - R)), proving (12111) once 
more. In fact, if this were not the case, we would have, like before, 

rRn 

j\E\ < R+Rn a r (dB r nE n ) dr = \E n n(B Rr , \B R+Rn ) \ < \E n nB R J-\E n nB R \ < l\E\ , 

which is a contradiction. This completes the proof of (12.31) . 

It remains to prove (j2.4[) (which implies (12.51) ) and, to do so, we again distinguish 
two cases. Assume hrst that E is bounded and note that 


F n AE\ — 2\F n \ E\ + \E\ — \F n 


(2.7) 
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Since F n D E u GBr, we have |F n | > \E n nBji\ and, by local convergence, \E h P\Br\ —> 
|E fl Br | = \E\. Thus, liminf \F n \ > \E\. On the other hand, we have F n C E n D B 2 R 
and therefore \F n \E\ < |( E n \E) D B 2 r\. By local convergence, | (E n \E) D B 2 r\ —> 0. 
Because of (12.71) we conclude that \F n AE\ —> 0. Finally, if K is compact, then 
\G n n K\ < \E n O (R‘ 1 \ Br ) fl K\. By local convergence, the latter converges to 
| E D (M d \ Br) fl K | = 0, so we conclude that G n —> 0 locally. 

Now assume that E is unbounded and note that 


\F n AE\ — 2| E \ F n | + \F n \ — \E\ . (2-8) 

Since F n C E n 0 Br u , we have \F n \ < \E n n BR n \ — \E\. On the other hand, we have 
\E \ F n \ = | (E \ E n ) 0 B rn | + \E \ B Vn \. By dominated convergence, since r n —> 00 , we 
have | E \ B r J\ —>■ 0. Now given e > 0 let p > 0 such that | E 0 (K d \ B p )\ < e. Then, 
if r n > p, \(E \ E n ) O B rn \ = \(E \ E n ) D B p \ + \(E \ E n ) O B r . n O (R d \B P )\. We have 
| (E \ E n ) HB p \ —)■ 0 by local convergence and | (E \ E n ) 0 B rn 0 (M d \ B p ) | < e. All this 
proves that | E \ F n \ —> 0. Because of (12.8p we conclude that \F n AE\ —> 0. Finally, if 
K is compact, then K fl (M d \ B rn ) = 0 for all sufficiently large n, so \G n DK\ =0 and 
we conclude that G n —>- 0 locally. This concludes the proof of the lemma. □ 


Clearly, the sets F n and G n from the previous proposition satisfy 

\F n \ + |G>i| = \E n \ . 

We next show that, at least asymptotically, some non-local functionals are also additive 
with respect to this decomposition. For 0 < A < d let 


UE) = - 


ExE 


dx dy 

\x ~ y\ x ' 


Lemma 2.3. Let F n and G n be sequences of measurable sets in W 1 with uniformly 
bounded measure and \F n D G n \ = 0 for all n. Assume that 


F n —$■ E (globally) and G n —> 0 locally 
for some set E. Then for all 0 < A < d, 

h(F n U G n ) = I x (F n ) + I x (G n ) + o(l). 

and 

I\(F n ) — I\(E) + o(l). 


(2.9) 


( 2 . 10 ) 


Proof. As a preliminary to the proof we note that for an arbitrary set F C M. d of finite 
measure, by a simple rearrangement inequality US Thm. 3.4], 


dy 


\x-y | 


a< 


dy 

w 


= C\F | 


(, d-X)/d 


with some explicit C depending only on d and A. Moreover, 

f dy 


\x-y\ 


->■ 0 


as \x —> 00 . 


( 2 . 11 ) 


( 2 . 12 ) 
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(This follows from [LL, Thm. 2.2] by decomposing |x|” A in a short and a long range 
part.) 

We now turn to the proof of (12.91) . We need to show that 


2/a(-Fti, G n ) 



dx dy 
x — y\ x 


-A 0. 


Using (12.lip and F n —)■ E we find 


\h(F n ,G n )-I x (E,G n )\ <I x (F n AE,G n ) < C\F n ^E\^ d \G n \ ^0, 


so it is enough to prove I\(E,G n ) —> 0. Let £ > 0 be given and use (12.12)1 to find 
R > 0 such that f E \x — y\~ x dy < e for |x| > R. We decompose and bound 

h(E, G n ) = Ix(E, G n n B r ) + I X (E, G n \ B R ) < C\E( d ~ x ^ d \G n fld R |+ e\G n \ B R \. 


Since G n —> 0 locally, we have \G n nB R \ —$■ 0 and therefore I\(E, G n ) —> 0, as claimed. 
The proof of (12.101) is similar. □ 


3. Application to the liquid drop model 


3.1. Results. In this section we return to the liquid drop model discussed in the 
introduction. Thus, we assume d — 3, A = 1 and abbreviate ii(fl) = D( fl). We 
recall that the binding energy of a set and the minimal binding energy with nucleon 
number A were defined as 


£(Q) = Per 12 + D(Q) and E(A) — inf{£(0) : |fl| = A} . 

It is easy to see (and proved, for instance, in [CO]) that 

E{A) < E(A') + E(A - A') for all 0 < A' < A. (3.1) 

Our first theorem states that if this inequality is strict, then there is a minimizer for 
E(A). 

Theorem 3.1. Let A > 0 such that 


E(A ) < E(A') + E(A - A') for all0<A'<A. (3.2) 

Then the infimum defining E(A) is attained. Moreover, any minimizing sequence has 
a subsequence which, after a translation, converges (globally) to a minimizer. 

We recall that (global) convergence of sets was defined at the beginning of Section [2] 
It means that the measure of the symmetric difference between the sets and their limit 
tends to zero. 

Of physical interest is the quantity 


the binding energy per particle. 



Theorem 3.2. There is an A* > 0 such that e(A*) = inf{e(A) : A > 0}. Moreover, 
the infimum defining E{A *) is attained. 
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In other words, there is a set hi* C M 3 with finite perimeter such that 


Per n + D(Q) ^ Per O* + T>(0*) 

[fij “ fO} 


for all C M 3 with finite perimeter. 


Theorem 13.21 has the following simple corollary. 


(3.3) 


Corollary 3.3. As A —» oo, e(A) e(A*). Moreover, e{kA *) = e(A*) for all k E N. 

Proof of Corolla, ry \3.?A Subadditivity ( 13.111 and non-negativity of E(A) imply by ab¬ 
stract principles that liixu-^oo e(d) exists. (This is sometimes called ‘Fekete’s lemma’.) 
Iterating ( 13 . 11 ) we infer that e(kA *) < e(A*) for all k E N. On the other hand, since 
A * is the global minimium of e(A), we have e(kA *) > e(A*). So e(kA *) = e(A*) for 
all k E N and, since we know that e(A) has a limit, this limit must be equal to e(A*), 
as claimed. □ 


Our final theorem is 


Theorem 3.4. The set {A > 0 : E(A) has a minimizer} is closed. 

This proves, in particular, that for A C2 = sup{d > 0 : E(A) has a minimizer} there 
is a minimizer. 


3.2. Proof of Theorems 13.1113.21 and 13.41 We now prove our three main theorems 
about the liquid drop model. 

Proof of Theorem \3.1[ Let (f2 n ) with |f2 n | = A be a minimizing sequence for E(A). 
Since Perf2 n < £(£l n ) = E(A) + o(l) is uniformly bounded, Proposition 12.11 yields 
a set 0 C M 3 with 0 < |0| < A and Perf2 < lim inf Per Q n such that, after passing 
to a subsequence and a translation, we have Q n —> 0 locally. Moreover, by Fatou’s 
lemma, we have D(Q) < lirninf D(Q n ), so 5(f2) < liminf£(O n ). Thus, fl will be a 
minimizer provided we can show that |Q| = A. Moreover, it is easy to see that Q n —> fl 
locally and |f2 n | —y \A\ implies that fl n —> globally. Therefore both statements of 
the theorem follow if we can prove that |Q| — A. 

We argue by contradiction and assume that |Q| < A. Then according to Lemma 12.21 
we can write D, n = F n U G n , F n D G n = 0, such that 

Per fl n > Per + Per G n + o(l). 

Moreover, F n —y globally and G n —y 0 locally which, by Lemma 12.31 implies that 

D(n n ) = D(n) + D(G n ) + o( 1 ). 


E(A) = £(n n ) + o{l) 

> £(Q) + £(G n ) + o( 1) 

>^(|n|) + ^?(|G n |) + o(i). 


Thus, 
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Since \G n \ = |12 n | — |F n | -A A — |12| and since A i —y E(A) is continuous on (0, oo) (in 
fact, A~ 2 / 3 E(A) = inf {Pern; + AD(u) : |o;| = 1} is concave as an infimum over affine 
linear functions), we obtain in the limit n —> oo 

E(A)>E(\n\) + E(A-\n\), 

which contradicts assumption 13.21 Thus, |12| = A. □ 


We prove Theorem 13.21 via the auxiliary minimization problem 

e<(A) := inf j^yp : 0 < |12| < aJ . 

The idea of relaxing the equality constraint |12| = A to |12| < A is reminiscent of |LS ], 
Lemma 3.5. For any A > 0, the infimum defining e<(A) is attained. 


Proof. Let (12 n ) with |12 n | < A be a minimizing sequence for e<(A). By the isoperi- 
metric inequality, we have e<(A) + o(l) = £{£l n )/\£l n \ > C\Ft n \~ 1 / 3 , which implies 
that liminf|12 n | > 0. Moreover, Perl2 n < £!(12 n ) < A(e<(A) + o(l)) is uniformly 
bounded, so as in the proof of Theorem 13.11 we obtain, after passing to subsequence 
and after a translation, a set 12 C R 3 with 0 < |12| < A and with Q n —y 12 locally and 
£’(12) < lim inf £’(12 n ). We now distinguish two cases, according to how 1121 compares 
to lim inf 112„ |. 

If |12| > lim inf 112 n |, then £1(12)/|12| < lim sup (£(Q n )/ 112„|) = e<(A), so 12 is a 
minimizer for e<(A). 

If |12| < lim inf |12 n |, then, as in the proof of Theorem 13.11 


e< (A) 


£(12) +£(G n ) 
1121 + | G n | 


+ o(l)> 


£ 1 ( 12 ) 

1121 + | G n 


+ 


\G n \ 

1121 + | G n 


e<(A) + o(l). 


Rearranging the terms, we obtain 1121 e<(A) > £1(12) + o(l), which again means that 
12 is a minimizer for e<{A). □ 


Proof of Theorem \3. 6 A Clearly, A i— y e<(A) is non-increasing and non-negative, and it 
is continuous because of the continuity of E i —> E(A) (see the proof of Theorem 13. II) . 
Let 

A * := sup {A > 0 : there is an A' > A with e<(A') < e<(AL)} . 

We claim that A* < oo. Clearly, this implies that e<(A) = e<(A*) for all A > A* and 
therefore that e(A) > e(A*) for all A > 0. Moreover, the minimizer 12* for e<(A*), 
which exists by Lemma 13.51 has |12*| = R* and therefore is also a minimizer for E(Af). 

We argue by contradiction and assume that A* = oo, that is, there is an increasing 
sequence (A n ) with A n — * oo such that e<(A n+ 1 ) < e<(A n ) for all n. According 
to Lemma 13.51 there are sets 12 n with <£(12 n ) = e<(A n )|12 n | and |12 n | < A n . The 
strict inequality e<(A n+ 1 ) < e<(A n ) implies that |12 n+ i| > A n and so, in particular, 
Ad n := |12 n | —> oo. Since the sets 12 n minimize e<(A' n ), they also minimize E(A' n ), but 
the existence of minimizers for this problem with arbitrarily large A' n contradicts the 
result in [LCH Thm. 2]. This proves that A* < oo, as claimed. □ 
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Remark 3.6. The minimization problem of Theorem 13.21 is equivalent to the following 
scale-invariant minimization problem, 

/ = inf | ^ Per --)^ D ( Q )' /3 . n cR 3 of finite perimeter 

l l“l 

(We know from [KM, Lem. 7.1] that I > 0.) In fact, in order to minimize £(12)/|12| 
we can minimize separately over shape and size of hi = £oj, that is, 


inf 

n 


£( 12 ) 


inf inf (i 1 Per a; + i 2 D(u)) 
M=i ^>o v ' 


For fixed u>, the inhmum is attained at = (Per uj/(2D(co)) 1 ^ 3 and we have 


inf = 2“ 2/3 • 3 • inf (Per uj) 2/3 D(u) 1/3 = 2~ 2/3 • 3/ . 

O |12| ]w|=l 

In particular, if uj has |co| = 1 and minimizes (Per 12) 2,/3 Zl(a;) 1//3 , then 12 := £ u cu 
minimizes £(12)/|12|. 


Proof of Theorem \3.f\ Let A n be a sequence with A n —> A e (0, oo) such that E(A n ) 
has a minimizer 12 n with |12 n | = A n . As observed in the proof of Theorem 13.11 
A t —y E(A) is continuous, so E(A n ) —y E(A) and Perl2 n is uniformly bounded. As in 
the proof of Theorem 13.11 we may assume that 12 n —> 12 locally for a set 12 C K 3 with 
0 < 1121 < A and £(12) < liminf £(12 n ), and it remains to prove that |12| = A. 

To prove this, we use a bound (implicitly) contained in [LOj . Namely, for every 
e > 0 there is a constant C e > 0 such that if A' > e and 12' C K 3 is a minimizer for 
E(A'), then 


diarn !2 / < C e A . (3.4) 

This is stated in |IKMI Lemma 7.2] with e — 1, but the same proof works for e < 1. 
(Technically speaking, in order to avoid ambiguities with sets of measure zero, 12' is 
replaced by the set {x G M 3 : limsup,.^ 0 112' D B r (x)\/\B r (x)\ >0}. By Lebesgue’s 
differentiation theorem, this set coincides with 12' up to sets of measure zero, so none 
of the terms in the minimization problem changes under this replacement.) 

Applying (13.41) to 12' = 12 n we infer that lim sup^^ diam 12 n < oo. Thus, the sets 
12 n are contained in a fixed ball, and then local convergence 12 n —y 12 implies global 
convergence, so 1121 = lim |12„| = lim A n = A, as claimed. □ 


Remark 3.7. Theorems 13.11 and 13.41 remain true, mutatis mutandis, for the functional 
£(12) = Per 12 + I a ( 12), 12 C K d , with d > 2 and 0 < a < d. (For the analogue of the 
diameter bound (13.41) see [ KM1 Lemma 7.2].) Theorem 13.21 and Corollary 13.31 remain 
true under the additional assumption 0 < a < 2. (Under this assumption the analogue 
of the Lu-Otto non-existence result remains valid, see [KMl Thm. 3.3].) 
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4. Further remarks about the liquid drop model 
4.1. A lower bound on A*. Let us consider the binding energy per particle for balls, 


,(ball) 


O) := - 4 ~ 1/3 1^273 + ■ 42/3 


Peril i0/ , D(B) 

|fi| 5 /3 ’ 


(4.1) 


where B is the unit ball in R 3 . (So \B\ = 47t/ 3, Peril = 47 t and D(B ) can be 
computed, but we don’t need this.) Clearly, e ( b “^(A) decreases up to some A^ alV> > 0 
and then increases. Setting the A-derivative of (A) equal to zero we find 

I ill Peril 


j^iball) _ 


(4.2) 


2 D(B) ' 

The main result of this section is the following quantitative bound on the curve e(A). 

Proposition 4.1. Let 

A 0 := sup {A 0 > 0 : Ak e(A) is strictly decreasing on (0, A 0 )} 

Then A 0 > A^ >all ' > with equality iff E(A 0 ) is minimized by a ball. 


Note that this gives, in particular, the lower bound A* > A 
from Theorem 13.21 


(ball) 


on the number A* 


Proof. We know from [KM] that E(A) is minimized by balls for small A and therefore 
e(A) is decreasing for small A, so Ao > 0. From Theorem 13.21 we know that e(A) has a 
global minimum, so Ao < oo. By assumption we have e(A) > e(Ao) for all 0 < A < Ao 
and, therefore, 

E(A') + E(A 0 - A') > A'e(Ao) + (A 0 - A , )e(A 0 ) = E(A 0 ) for all 0 < A' < A 0 . 

According to Theorem 13.II this implies that there is a minimizer Oo for -F(Ao). 

The next step is to derive a virial relation for Oo- Since e(A) has a local minimum 
at Ao, the function 

£(«U) .- Terfio f 2 -P(W) 

I COo | | Oo | | Oo | 

has a local minimum at i = 1. Setting the derivative at i = 1 equal to zero we 
conclude that 

Perflo = 2D(fl 0 ). (4.3) 

We now use (14.31) to prove a lower bound on Ao- For the proof we will use the 
inequalities 

Per E > Per E* and D(E*) > D(E ), (4.4) 

where E* denotes a ball of the same measure as E. The first inequality is just a 
rewriting of the isoperimetric inequality and the second one follows, for instance, from 
the Riesz rearrangement inequality [LL , Thm. 3.7]. From these inequalities and (14.31) 
we deduce that 

ElE d(SIo) ^ , 


1 Ppr R 

D(n 0 ) = -Perfio >-PerfR =-—— D(Q*) 2/5 > -- , NO/K . 

v o; 2 u ~ 2 0 2 D(B) 2 E v oJ ~2 D(B) 2 E 
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(The middle equality here just uses scaling.) Thus, 


D(fin 1 3/5 > - 


1 Per B 


0) ~ 2 D(B) 2 / 5 

On the other hand, again by (14.41) . 

D(n 0 ) 3/5 < d(q *) 3/5 = ®nsi = D( ,^, ) | 3/5 |Ool • 


IBI 


(The middle equality again just uses scaling.) Combining the last two inequalities and 
recalling (j4.2|) we obtain 


Ifiol > 


B | Per B 
2 D(B) 




as claimed. 

Finally assume that f2 0 is not a ball. Then the strict rearrangement inequality 
D(Sl‘ 0 ) > D(Qq) [Llj implies, by the same argument as before, the strict inequality 
| O 0 1 > A^ all \ This concludes the proof. □ 


Remark 4.2. The same proof shows that the virial relation (14. 3 j) holds whenever Aq is 
a local minimum of e(A) and E(Aq) is attained by some Oq- 


4.2. Dissociation into balls. Let us consider 

e(A) := A~ l inf {£(£)) : |0| = A , is a countable union of disjoint balls} . 

This is the same as the energy e(A) per particle except that we restrict the allowed 
sets to be countable disjoint unions of balls. Clearly, it suffices to consider the case 
when the individual balls are infinitely far apart and therefore 

f oo . 00 

e{A) = inf : ^A k = A 

l. k =1 k =1 

with the energy e^ ball \Ak) from (14.11) . It turns out that this minimization problem can 
be almost explicitly solved. 


Proposition 4.3. For any A > 0, 

e(A) = min e^ ball) (A/k). 
fce N 

That is, e(A) coincides with the energy per particle of K infinitely separated balls 
of volume A/K, where K is to be optimized over. The energy per particle of such 
a ball satisfies e^ ball \A/K) ~ [\dB\/{A l ^\B\ 2 ^))K 1 ^ as K — > oo and therefore the 
inffinum over K is indeed attained at some finite K. With a little more work one 
can compute explicit numbers 0 = a 0 < cq < a 2 < a 3 < ... with —* oo such that 
e(A) = e ( 6a C (A/k) for ak-\ < A < ak- For instance, 

_ 2 — 2 2 / 3 \B\ Per B 
ai ~ 2 2 /3 - 1 D(B) ' 

The proof of Proposition 14.31 is based on the following technical lemma. 
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Lemma 4.4. Let A > 0 and f{9 ) = (6 2 / 3 + (1 — 6 1 ) 2 / 3 ) + A(# 5 / 3 + (1 — 6 ) 5 / 3 ). Then for 
all 6 G [0,1] 

f(0) > min{/(0),/(l/2),/(l)} 
with strict inequality unless 6 G {0,1/2,1}. 


Proof of Lemma \4-4\ Differentiating we find 

- S(«)) , (4.5) 

where 

= (1 - Op 3 - 9 2 / 3 

' 1 ' 9-1/3 _ (1 _ (lj-l/3 ' 

Clearly, g is symmetric with respect to 9 = 1/2 and satisfies g(0) = p(l) = 0. By 
straightforward differentiation we find that g is strictly concave, so it attains its max¬ 
imum at 6 = 1/2 and we compute g( 1/2) = 1. Thus for every 0 < 7 < 1 there is a 
unique 9( 7 ) G [0,1/2) such that g{9) = 7 iff 0 = ± 0 ( 7 ). 

Returning with this information to (14. 5ft we find that, if 2/(5A) > 1, then /' has a 
single zero at 6 = 1/2 and, if 2 /( 57 ) < C then /' has three zeros 6 = 1/2, ±0(2/(5A)). 
Since / is clearly increasing near 0 and decreasing near 1, the single zero of f at 
6 = 1/2 for 2/(5A) > 1 and the two zeros at 9 = ±9 (2/(5A)) for 2/(5A) < 1 must 
correspond to local maxima of /. Thus, / attains its minimum at 6 G {0,1} for 
2/(5A) > 1. For 2/(5A) < 1, / has a local minimum at 9 — 1/2 and it attains its 
minimum on a subset of {0,1/2,1}. This proves the lemma. □ 



Proof of Proposition \ 4 ■ -4 We only need to prove the inequality >. We first note that 
in the definition of e(A) we can restrict ourselves to finite sums. Indeed, for any ARN 
put A k = Y,T=i<Ak, so that 


A 


^ A 

k =1 


^ e^ ball \A k ) > 



3 ( ball) 


(A k ) + ^ ball \A K ) 


^ e^ ball \A 


k 


Since ae ( ' ball \a) —> 0 as a —* 0, the last term goes to zero as K -A 00 . The first term 
can be bounded from below by the claimed expression if we prove the result for finite 
sums, as we will do now. 

In fact, we will prove that 


K 


Y. A ^ = A 


nun e 

K>k&N 


(ball) 


(A/k) 


(4.6) 


. fc=l 


k= 1 


by induction over K e N. For K = 1 is clearly true. Now assume that K > 2 
and that the assertion is proved for 1,..., K — 1. The set {(All,..., Ak) G [0, A] K : 
Yf, = A} is compact and the function to be minimized on the left side of (14.6p is 
continuous (note again that ae^ ball \a) is continuous at a = 0), so there is a minimizer 
Ah°h If A ^ = 0 for some k, then (14.61) follows by induction assumption. Therefore we 
may assume that A ^ ^ 0 for all k. Let 1 < i < j < K be arbitrary. We shall show 
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that Af ] = A^\ which will prove (14.61) . The contribution from these two numbers to 
the left side in (14. 6 p is 


A 


(o) 


A 


e ' Mi >(A< 0) ) + = (e 2/3 + (i - ») 2/3 ) (4°> + .4f ) 2/s AIAL 

+ (e 5/3 + (1 - ff) s/3 ) (A 0) + ^ 0) ) 5/3 


with 9 = A^/(A^ + A^). By construction we have 6 7 ^ 0,1. If we had 9 7 ^ 1/2, 
then by Lemma [4.41 this contribution would be strictly less than for 9 = 1/2. Thus, 
replacing A^ and A^ both by (A® + A < J °' l )/2 we would get a strictly smaller energy 
while preserving the constraint, but this contradicts the minimality of A*- 0 / Thus, we 
have 9 = 1/2, that is, A = A^\ which is what we wanted to prove. □ 
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